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Abstract. - 

The quantum-limited linewidth of a laser cavity is enhanced above the Schawlow-Townes value 
by the Petermann factor K, due to the non-orthogonality of the cavity modes. The average 
Petermann factor (K) in an ensemble of cavities with chaotic scattering and broken time-reversal 
symmetry is calculated non-perturbatively using random-matrix theory and the supersymmetry 
technique, as a function of the decay rate V of the lasing mode and the number of scattering 
channels N. We find for N 3> 1 that for typical values of V the average Petermann factor 
(K) oc vTV 3> 1 is parametrically larger than unity. 



The study of resonant scattering goes back to the early work of Breit and Wigner |lj] , and 
was developed extensively in the context of nuclear physics [|| . The Breit- Wigner resonance 
is described by a frequency-dependent scattering matrix S{uj) with elements 

Sum = Snm + &n a 'm{u ~ ^ + iT/2)~ 1 , (1) 

where ft and T are, respectively, the center and the width of the resonance, and a n , a' m are the 
complex coupling constants of the resonance to the scattering channels n, m. (In the presence 
of time-reversal symmetry the scattering matrix is symmetric, hence o~ n = a' n .) The coupling 
constants of a Breit- Wigner resonance are related to its width by the sum r ule @ 

5>XJ 2 = r 2 . (2) 

The physics behind this sum rule is that narrow resonances are quasi-bound states requiring 
a weak coupling to the continuum. 

The Breit- Wigner theory holds in the limit T — > only. More generally, one may consider 
a pole Q — iT/2 somewhere in the lower half of the complex w-plane and ask for the form of 
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the scattering matrix near that pole. The separable form (|l|) remains valid, even if T is not 
small. The sum rule however, should be replaced by 

5>X„l 2 = ^r 2 , K >l ( 3 ) 

The total coupling strength of the resonance is larger than required by its width by a factor 
K. This factor has been studied theoretically and experimentally in laser cavities |5|, q, [|, 
where it governs the quantum-limited lincwidth, 

Suj = K^, (4) 

of the radiation from the lasing mode (with output power /) . The lincwidth is larger than the 
Schawlow-Townes value T 2 /2I || by the factor K. In this context K is called the Petermann 
factor or excess noise factor (E). 

The optical cavities considered in these studies typically have a simple shape, that is 
characterized by integrable rather than chaotic wave dynamics p0|| . Recently, three of us 
investigated the statistics of the Petermann factor in the presence of chaotic scattering j|. 
That work was restricted to the case of a single scattering channel (N = 1), meaning a scalar 
scattering matrix. If the cavity is coupled to the outside world via an opening of area A, 
then a single scattering channel requires A < ^ 2 /27r, with A the wavelength of the radiation. 
The factor K is only slightly larger than 1 for such an almost closed cavity, and measuring 
the deviation from unity would be quite problematic. For a larger K one has to increase A, 
thereby increasing the number of scattering channels N ~ 2irA/\ 2 . 

The fundamental question addressed in this paper is how K increases with increasing N . 
We present for the case of broken time-reversal symmetry a non-perturbative calculation of 
the average Petermann factor (K) in an ensemble of chaotic cavities, valid for arbitrary N and 
T. A non-perturbative calculation is necessary because the mode selected for laser action has a 
decay rate T of the order of the typical smallest rate To = TNA/2ir, with A the mean spacing 
in frequency of the resonances and T the transmission probability of the channels through the 
opening. (The time 1/To is the mean dwell time of photons in the cavity.) Existing methods 
of large- N perturbation theory require that V > Tq, and are therefore not applicable to 
this problem. Our non-perturbative method is based on the supersymmetry technique jl^, [l3| , 
applied to eigenvector correlations of non-Hcrmitian random matrices. Our conclusion is that 
(K) oc y/~N ^> 1 can become parametrically large for N 3> 1. The square-root increase is the 
outcome of a subtle balance between the smallness of T and the largeness of N, and could not 
have been anticipated from perturbative arguments, which would imply a linear increase. 

To formulate the problem within the framework of random-matrix theory we proceed as 
in ref. Q. The modes of a closed chaotic cavity, in the absence of time-reversal symmetry, 
are eigenvalues of an M X M Hcrmitian matrix H with independent Gaussian distributed 
elements [|l0[ [l4| . The limit M — > oo at fixed spacing A of the modes is taken at the end of the 
calculation. The opening in the cavity is described by an M x N coupling matrix W. Modes 
of the open cavity are complex eigenvalues of the non-Hermitian matrix H = H — iirWW^ 
which determines the scattering matrix ||, 

S = 1 - 2mW\uj - H^W . (5) 

The matrix H, — [/diag (z l5 z 2 , . . . , zm)U~ 1 has complex eigenvalues Zj and the matrix U 
contains the right eigenvectors of 7i in its columns. (The matrix U" 1 contains the left 
eigenvectors.) Because Ti is not Hermitian, U is not unitary and the left or right eigenvectors 
are not orthogonal among themselves. If the cavity is filled with a homogeneous amplifying 
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medium (amplification rate l/r ) one has to add a term i/2r a to the eigenvalues, shifting them 
upwards towards the real axis, while the eigenvectors are unchanged. The lasing mode is the 
eigenvalue zi = fl — iY/2 closest to the real axis, and the laser threshold is reached when the 
decay rate T of this mode equals the amplification rate l/r a . Near the laser threshold it is 
only necessary to retain the contribution from the lasing mode. Combining eqs. (|l|), (^), and 
(||) one finds that the Petermann factor is given by [[|, [| 

K = (C/t U) u (U- 1 U-^) u . (6) 

This formula shows that deviations of K from unity result from the non-unitarity of U , hence 
from the non-orthogonality of the cavity modes || . 

We will calculate the average Petermann factor at a given value of T and Vt. Following refs. 
[ fill [[U > we write this conditional average in the form 

(K) = (j2^ U ^n(U- 1 U-^) n 6(z-z 1 )){Y^6(z-z J ) 

3 
1 



J 1 

2\ 



- irp(z)H\{ tT (z-H)(z*-Hr)+ei) / < (?) 

where z = O — iT/2. The eigenvalue density p(z) — (J^ . 8{z — zj)) has been determined by 
Fyodorov and Sommers |Tj| , 

p(z) = (2vr/A 2 )^ 1 (^r/A) ^ 2 (^r/A) , (8) 
The parameter g > 1 is related to the transmission probability T of a scattering channel by 

.</ / i Q 

The evaluation of the quantity ({?]) proceeds by the supersymmetry approach. We introduce 
the generating functional (det(A — J)/ det(A + J)), where J is a source matrix and 

A=( A A n A A 12 ) = ( * t Z ~ H ) ■ (10) 
V A 2i A 22 J \ z* -W is J y 1 

The quantity np(z) (if) = — (tr[(A~ 1 )n] tr[(A _1 ) 2 2]) can be expressed in terms of derivatives of 
the generating functional with respect to J. Following the technique of supersymmetry |TJ, 13 
one can express the Petermann factor as an integral over a 4 x 4 matrix Q, 

p{z)(K) = -JL lim [ dQe- c ^ Str {Qa x P P+) Str {Qa x P P-) , (11) 

£(<2) = Str (Qcr x ) - — Str (a z Q) + NStr ln(l + wQa z ) , (12) 

where w is related to T by T = 4w(l + w)~ 2 . The matrix Q obeys the non-linear con- 
straint Q 2 = 1 and belongs to the coset space of the unitary non- linear er-model |l2j. It 
is a "supermatrix" , meaning that it has an equal number of commuting and anticommuting 
variables. The symbol "Str" denotes the graded trace of a supermatrix, Str A = tr PqA, where 

( 1 ) We assume for simplicity that each channel has the sa me q . Expressions for for channel- 
dependent g n can be found in ref. jl3|], and our main result (|lj) given below remains valid for the 
general case. 
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Pq = diag(l, — 1, 1, — 1). The Pauli-matrices a Xl a z are embedded in the 4x4 supermatrix 
space as a z = diagfl, 1, — 1, -1) and {a x )jk = 5j,fc+a + &j+2,k- Furthermore, P± = |(1 ± a z ). 

The integral (11) can be evaluated by int rodu cing Efetov's parameterization of the matrix 
Q for the unitary er-modcl |L3]. The result it ( 2 ) 



p(z)(K) 



7T d 

■— r hm e— 

2A 2 e^o de 



Ai + A 2 ( g + A 2 
5 + Ai 



dA2 > i 

1 Ai — A2 



N 



X Jo 



A 2 — 1 I exp 



wT (X 



A 2 ) 



(13) 



where Jo is a Bcssel function. The remaining two ordinary integrals can be carried out as in 
the corresponding calculation of p(z) jl3|. We finally obtain 



(K) 



1 



2S(y) 



S(y) = - [ dx T x (x) ^- T 2 (x) , y = ttT / A 
Jo dx 



A more explicit expression can be obtained by successive partial integrations: 



S(y) 



N-l . 

E- 

k=0 



1) 



■ V 



-gy _d_ sinhy 

dy y 



(14) 



(15) 



(if) = 1 + 27rr/3gA. 



For N = 1 and r « A we recover the single-channel result of ref. 
Eq. ( |l4[) generalizes this result to any value of N and T. 

To study the limit N 3> 1, it is convenient to derive integral expressions by replacing 



sinh y/y in Ti with i j -^dXt 
the result 



Ay . This was done in ref. fl3]] for the eigenvalue density, with 



2y*(N-l)\ i M 



(9+1)9 



dtt N e~ l 



(16) 



Similarly, we find for the function S'(y) that determines the Petermann factor the expression 



i 



4y 2 (JV-l)lVi) 



d< [*-07-l)»][t-(ff + l)tf]t 



N-i e -t 



(17) 



For simplicity we restrict ourselves in what follows to the case T = g = 1 of fully transmitted 



scattering channels. When N ^> 1, the factor t 



JV 



has a sharp maximum at t = N and one 



can apply the saddle-point method to evaluate the t-integration. This gives the eigenvalue 
density || p(z) w iV/7rr 2 for T > T (with T = NA/2ir for T = 1). The density vanishes 
exponentially as T drops below r . Application of the same saddle-point method to eq. JTe| ) 
gives (K) w (27r/A)(r-r ). If T-To = O(r ), then this estimate would imply that (K) cx N. 
This linear scaling with N breaks down in the cutoff regime T < To, which is precisely the 
relevant regime for the laser cavities. 

In order to study the cutoff regime we define the rescaled decay rate u = \J 'N /2(T /Tp — 1) 
and take the limit N 3> 1 at fixed u. Expanding the integrands in eqs. ([t6]) and (17) around 
the saddle-point and keeping the first non-Gaussian correction, we obtain 



(K) = V2N[F{u) + u}+ F(u) [2u + § u 3 + |(1 + u 2 ) F{u)] + 0(N~ 1/2 ) 



(18) 



( 2 ) For r real positive the integrals (O) and ([TJ|) are formally divergent. The proper procedure for 
the further evaluation is to choose T imaginary and to perform the analytic continuation at the 
end of the calculation. 
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where we have abbreviated F(u) = 7r -1 / 2 exp[— u 2 } [1 + erf(it)] -1 . The eigenvalue density in 
the cutoff regime has the form p = (27r/7VA 2 )[l + erf (u)} + 0(A^ 3/2 )- For u = (hence 
T = To) eq. ( |l8|) simplifies to 

l2N 4 



while for u < -1 (hence V - T < -ViVA) we find 

(^-i = ~Jf = A. (20) 
u V 2 1 — 1 

To determine how (-ftf) scales with N we need to estimate the most probable value of T 
for the lasing mode. The decay rate Ti of the lasing mode is the smallest among the cavity 
modes that are amplified by stimulated emission. The number m a of amplified modes is 



3> N in the "good cavity" regime Q. In a simplified model |L6|, we can determine T; as 
the smallest of m a independent random variables Ti, . . . ,T ma with probability distribution 
P(T) = (A/2)p(z) S3 (l/2r )[l + erf («)]. The distribution of f ; is then given by 



Pi(T l )=m a P(Ti) 



1 - \ dTP(T) 
Jo 



This distribution has a pronounced maximum at a value u max determined by 



[1 + erf(u max )] 2 2 V 2iV 



(21) 



(22) 



If m a is comparable to vN, we have w m ax = 0(1) and the Petermann factor is to close to its 
value at u = 0, eq. ( | l9]). For m a y/N (which is realized in the good-cavity regime) we find 

"max ~ -\Jlii(m a /VN) ~ -V m m a < -1 and 



(23) 

This is still parametrically larger than unity as long as m a <C e N . 

We now compare our analytical findings with the results of numerical simulations. We 
generated a large number of random matrices Ti with dimension M — 120 (A/ = 200) for 
N = 2,4,6,8 (N = 10, 12). Figure | shows the average of K at given T. We find excellent 
agreement with our analytical result (|l4|). The scaling (K) oc y/N at T = To is shown in 
fig. 0. The inset depicts the distribution of at T — Tq for = 10, which can only be 
accessed numerically. We see that the mean Petermann factor is somewhat larger than the 
most probable value. 

So far we only discussed the case of broken time-reversal symmetry, because this greatly 
simplifies the analysis. From a numerical simulation with real symmetric H (shown as well in 
fig. ^) we find that the average Petermann factor for T = To is parametrically large also in 
the presence of time-reversal symmetry. The average increases again sublinearly with JV, but 
with an exponent that is smaller than in the case of broken time-reversal symmetry. 



( 3 ) A "good cavity" has a smallest typical decay rate (inverse mean dwell time) To that is small 
compared to the amplification bandwidth Q a . Since m a — fi a /A and To = NA/2n, it follows that 
m a > N. 
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In summary, we have calculated the average Petermann factor (K) for a chaotic cavity in 
the case of broken time-reversal symmetry. We find that for typical values of the decay rate 
r the average (K) cx y/N scales with the square root of the number of scattering channels N, 
hence with the square root of the area A of the opening of the cavity. This result could only 
be obtained by a non-perturbative calculation because laser action selects a mode with an 
untypically small decay rate. From a numerical simulation we find that the sublinear increase 
holds also in the case of preserved time-reversal symmetry. The quantity K deserves interest 
not only because it determines the quantum limit of the linewidth of laser radiation, via eq. 

. It is of fundamental significance in the general theory of scattering resonances, where it 
enters the width-to-height relation of resonance peaks and determines the scattering strength 
of a quasi-bound state with given decay rate T, via eq. (^|). Because K \s a, measure for the 
non-orthogonality of the eigenmodes, this provides a way to obtain additional information 
about an externally probed open system. 

This work was supported by the Dutch Science Foundation NWO /FOM and by the TMR 
program of the European Commission. The authors thank the Max Planck Institute for Physics 
of Complex Systems in Dresden for its hospitality during the workshop: Dynamics of Complex 
Systems at which part of this work was done. 
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Fig. 2. - Average Petermann factor at V = To as a function of N. The exact result Jl4[ ) for broken 
time-reversal symmetry (solid curve) is compared with the result of a numerical simulation (open 
circles). Also shown is the simplified expression (113) for large N (dashed curve). The solid circles are 
the numerical result for preserved time-reversal symmetry. The inset shows the distribution of K at 
r = To for N = 10 and broken time-reversal symmetry, following from the numerical simulation. 
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